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Abstract 


Thermobaric explosives (TBXs) have been primarily used for 
their blast, rather than for their fragmentation, characteristics. 
This work reports on an investigation, using a flash X-ray imaging 
technique, of the ability of TBXs to shatter metal casings and to 
propel the resulting fragments. Three casing materials were used, 
AISI 1026 steel, ductile iron, and grey cast iron, while two 
different TBX compositions were used, with C4 serving as a 
benchmark. The fracture behavior of the casings, as a function of 
explosive fill and material characteristics, was as expected. One 
TBX formulation exhibited a run distance to detonation. The 
Gurney equation was used to correlate and compare the final 
fragment velocities. It was found that a larger fraction of the 
explosive energy was available to propel fragments in these two 
TBX compositions than a comparable amount of C4. This fraction 
of energy was influenced by the confinement of the detonation 
products and the ignition delay of the metal powders. These two 
factors had a greater influence on the fragment velocities than did 
material characteristics. 


Keywords: Thermobaric Explosives, Fragment Velocity, Materials, 
Gurney Equation 


Introduction 


Aside from chemical weapons, warheads are designed to 
create either fragments or blast shock waves as their primary 
damage mechanism. Thermobaric explosives (TBXs) have 
been predominantly used in a blast role, due to their 
enhanced blast effect, which is a direct result of the 
secondary combustion of additives. The shock wave gen- 
erated by the detonation of TBXs is of a lower amplitude, 
but longer period, than that of conventional secondary 
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explosives. The advancing shock wave is sustained as the 
additives ignite and burn after the main detonation reaction. 
This sustained shock wave makes TBX detonations very 
effective in enclosed spaces. 

The release of energy from the detonation of conventional 
secondary explosives is essentially instantaneous and results 
in a very rapid increase of temperature and pressure. TBXs 
combine conventional secondary explosives with a high 
volume fraction of additives, such as aluminum. Because the 
metal additives are considered inert in the reaction zone, 
due to their protective oxide films, and the additives replace 
some of the explosive volume, less energy will be produced 
in the TBX’s primary detonation reaction [1]. Consequently, 
the initial temperature and pressure produced are lower for 
TBXs than for equivalent conventional secondary explosive 
charges. 

In previous work with three TBX compositions [2], it was 
found that the fragment mass distributions and fragment 
velocities could be modeled as conventional explosives of 
similar detonation velocity (VOD) and density, despite 
evidence of additive combustion (oxidation) within the 
casing. Most of the energy from the combustion of the metal 
additives used was released when the additive particles 
came into contact with atmospheric oxygen. This paper 
reports on an investigation, using a flash X-ray imaging 
technique, of the ability of two TBX compositions to shatter 
metal casings and to propel the resulting fragments. The 
casing dimensions are based on the PMR-2A landmine. The 
results are indicative of the fragment velocities attainable 
from a dual-purpose (blast/fragmentation) warhead. 

Fragment velocity prediction originated with the work of 
Ronald W. Gurney in 1943 [3]. Gurney developed some 
simple and effective relationships that linked the masses of 
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casings and explosives to the observed fragment velocities. 
The key assumption in Gurney’s derivation was that the 


Table 2. Test matrix indicating the cylinder materials and wall 
thicknesses. An “X” indicates a series tested. 


fraction of energy of the explosive used to propel the Explosive AISI 1026 steel DCI GCI 
fragments was the same for all situations of explosives and 6mm 6mm 
: : ; 3.8mm 6mm 
casing materials. Gurney determined that for a fragment 
released from a cylindrical shell: C4 xX 
TBX 1 xX xX x x 
a TBX 3 xX 
M 1\? 

V=V2E (3 + 5) (1) 


where V is the fragment velocity in m/s, M is the mass of the 
casing, C is the mass of the explosive (in the same units as 
M), and V2E is referred to as the Gurney constant and is 
quoted in ms '. Gurney’s work has since been adapted to fit 
many different casing geometries and to account for other 
detonation phenomena [4]. However, the simple relation- 
ship shown in Equation (1) was appropriate for the 
estimation of the maximum fragment velocities of the 
casings studied in this research. The Gurney constant is 
dependent upon the explosive and casing material, and has 
been documented for explosives in mild steel [5]. 


2 Experimental 


The pertinent data for the two thermobaric compositions 
used in this work are listed in Table 1. The monopropellant 
was a liquid hydrocarbon fuel with modest detonation 
properties. 

The values for the detonation velocities for these TBX 
compositions were determined experimentally. The values 
for /2E are approximated on the basis of a correlation 
similar to those found in Cooper [5]. For TBX 1 in a brittle 
grey cast iron cylinder, a prompt Gurney constant is 
required and estimated on the basis of data tabulated in 
[6]. The prompt Gurney constant is 1.53 km s"! (further 
details are in the appendix). 

The test matrix to determine the fragment velocities is 
presented in Table 2. For this series, each test consisted of 
two explosive trials. The X-radiographic (X-ray) images 
were taken to estimate the maximum fragment velocities 
and qualitatively compare the casing break-up of several 
combinations of explosive, wall thickness and casing mate- 
rial. 

All the cylindrical casings had the same internal diameter 
of 44 mm, a length of 102 mm, and a wall thickness of either 
3.8 or 6.0mm. A scale drawing is shown in Figure 1. The 
materials were steel tubing, ductile cast iron (DCI) and grey 


Table 1. Properties of the explosives used in this work. 


Table 3. Typical mechanical values of AISI 1026 steel, ductile cast 
iron and grey cast iron (Askeland and Phulé, 2003) [7]. 


Metal Ultimate Yield Elongation 
tensile strength 
strength 
AISI 1026 steel 490 MPa 410 MPa 15% 
DCI (65-45-12) 450 MPa 310 MPa 12% minimum 
GCI (Class 40) 280 MPa - 0% 


cast iron (GCI). These irons were included in order to 
ascertain if casing ductility affected fragment velocity. A 
summary of the mechanical properties of these materials [7] 
can be found in Table 3. 

The end caps were made to match the material of each 
cylinder. The top caps were made separately and brazed in 
place. A polypropylene test tube, less than 1 mm thick, was 
epoxied in place, once the top cap had cooled from being 
brazed onto the cylinder. This ensured that the liquid TBX 
fills would not leak out of the top cap while the casing was 
inverted during the loading of explosives. It also reserved 
the required volume inside the casing for the booster charge 
that would be placed into the test tube. 

To prepare the TBX 1 charge, granular magnesium was 
added to a casing until almost full. The monopropellant was 
then slowly added until saturation. The casing was tapped, in 
order to release any air bubbles, and excess monopropellant 
was removed from the surface. The casing was weighed at all 
stages to confirm the quantities of magnesium and mono- 
propellant. For the preparation of TBX 3, the physical 
mixing of ingredients was done by machine in an evacuated 
area. A mass of 120 g of this explosive was added to each 
cylinder. Once filled, the base cap was then attached using 
an epoxy adhesive. When the adhesive was sufficiently 
cured, the charge was transported to the trial area where the 
C4 booster and detonator were added. 

The booster charge of 20 g of C4 was inserted into the test 
tube and tightly packed. The RP-83 exploding bridgewire 


Name Composition Detonation Velocity Density Gurney Constant 
km s! gcm> ms! 

C4 RDX/plasticizer (91/9) 8.09 1.77 2.83 

TBX 1 monopropellant/magnesium 4.6 1.18 1.75 

TBX 3 monopropellant/aluminum/RDX 5.9 13 2.11 
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RP-83 Detonator 


C4 Booster 


C4 or TBX Explosive 


Figure 1. Scale drawing of a 6mm casing. Note that for these 
trials, the cylinder is placed horizontally, with the detonator to the 
right. 


detonator was directly inserted into the C4, and a small 
amount of inert putty was added to plug the remaining few 
millimeters of the test tube. 

The X-radiographic (X-ray) images were taken at very 
short intervals after the detonation. The test site is 
illustrated in Figure 2 and pictured in Figure 3. As can be 
seen in these figures, the charge was placed horizontally ona 
polystyrene riser. Steel blocks were placed to protect the 
X-ray film cassette and surroundings, and to ensure that the 
two radiation sources did not overlap on the film. 

Two X-ray images were taken at each trial, where the two 
X-ray sources were set for different delays. As each test had 
two identical trials, four images were collected for each 
explosive-casing combination. The images were then used to 
calculate the velocity achieved by the fragments. This 


Film cassette 


Metal basin 


X-ray source | 


X-ray source 2 


Charge 


Plywood charge box 


Figure 2. Flash X-ray site apparatus, courtesy of Tom Storrie, 
DRDC-Suffield. 
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Figure 3. Photo of Flash X-ray site and (a) lead wires. Examples 
of damage to blocks from previous firings can be seen as (b) 
casing wall impact, and (c) casing cap impact. 


velocity was then compared to that calculated by the 
Gurney equation. 

There are many ways in which the time-distance informa- 
tion from the images can be interpreted. The simplest is to 
take the slope of the line from a plot of distance versus time, 
yielding an average velocity. Another way is to assume an 
exponential relationship for the velocity with respect to time 
of the general form: 


V=Vy (1 7 exp(—- (3) (2) 


where V, is an asymptotic velocity reached after an infinite 
time t. Note that at t=t, the velocity is 0.632 of the 
asymptotic velocity. The parameters n and t are constants 
fitted to the experimental data, t has units of time. To use 
this equation directly, the value of n must be specified; 
commonly n= 1 [8]. The equation is integrated starting at 
time t= a and V= Oto give the distance travelled, d, at time ¢: 


t 
d=Vi(t4+re'e—a— re“) (3a) 
For a value of n = 2, the integrated equation is then: 
1 t 1 
d=V, t—5t aert (4,) —a+5t nerf (4/,) (3b) 


where erf is the error function, and a, t, and V, are 
determined by a least squares fit to the experimental data. 
Although an equation with three parameters is fit to four 
data points, the parameters can be interpreted as having a 
physical meaning. The parameter a is the time of the start of 
casing expansion and Tt is a time constant related to the 
quickness of acceleration. For the analysis that follows, the 
focus will be on Equation (3a), where n= 1. 


3 Results 
3.1 Run Distance 


X-ray images showed that TBX 1 did decompose in a 
detonation reaction as the casing expanded uniformly. 
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Figure 4. Expansion of a 6mm AISI 1026 casing filled with 
TBX 1, 80 us after initiation from the right. 


Figure 5. Expansion of a 6mm AISI 1026 casing filled with 
TBX 3, 60 us after initiation from the right. 


Figure 4 is a typical image for this explosive fill, here in a case 
of AISI 1026 steel with 6 mm wall thickness. However, the 
image of the TBX 3 explosive in the 6 mm AISI 1026 steel 
casing showed asymmetric expansion at 60 microseconds 
after initiation (see Figure 5). This behavior is interpreted as 
the consequence of the run distance, where the shock wave 
from the booster must travel a distance comparable to the 
radius of the casing before the explosive achieves steady state 
detonation. It may also be indicative of a critical diameter 
effect for this composition and geometry, although this is 
doubtful, as TBX 1 has a lower velocity of detonation, and it 
detonated uniformly. Thus a run distance of about the 
cylinder radius (20 mm) is required for TBX 3 when initiated 
with a C4 booster. For TBX 1, the run distance is less. The run 
distance significantly affected the mass distribution results 
from the TBX 3 explosive, as will be discussed in a following 
paper. There was no evidence of spalling for any of the 
explosives or materials used in these trials. 


3.2 Fragmentation Characteristics 


The X-ray images also confirmed much of the expected 
casing break-up behavior. C4 initially generated more well- 


defined and tightly packed axial cracks than TBX 1. The C4- 
filled casings also expanded more rapidly and uniformly 
than those casings containing either TBX fill. For a TBX 1 
fill, changing the wall thickness from 6mm to 3.8mm 
increased the number and definition of the axial cracks 
generated, and also led toa more rapid casing expansion and 
smaller fragments. Casings of AISI 1026 steel initially 
showed predominant axial crack formation, as seen in 
Figure 4. The DCI casings showed these axial cracks as well; 
however, there were cracks in other directions. The trend 
away from uniform axial cracking increased when consid- 
ering the GCI casing, the most brittle material tested. There 
is anoticeable difference in average fragment size and shape 
as a function of casing material. In all, the fragmentation 
trends are as expected. 


3.2 Fragmentation Velocity 


Table 4 presents the calculated velocities from Gurney’s 
equation and the results of the two models above. The 
calculated velocities from the Gurney equation have an 
assumed error of 12%. The fit of the data to the linear 
equation, as seen in Figures 6 through 8, is good, indicating 
that the expansion of the cylinder occurs at nearly constant 
velocity. The exponential Equation (3a) form = 1 resulted in 
higher final velocities in all cases, and the largest discrep- 
ancy between these two equations is with TBX 1 in the DCI 
and AISI 1026 casings with a 6 mm wall thickness. It should 
be noted that the data were also fit to the exponential 
Equation (3b) with n=2. In this case, there was no 
appreciable difference in the final velocities determined 
using this equation and the linear equation. There was also 
no appreciable difference in the value of the time constant, 
t, forn = 1 or for n = 2, with the exception again of TBX 1 in 
the two casings noted above. The benefit of the exponential 
fit is, then, in the value of the time constant, giving an 
indication of the acceleration period of the fragments. 

The sum (a+) indicates the time during the trial where 
the velocity of the casing is 0.632 of the final velocity. 
Considering the results with C4 in Table 5, this sum is about 
40 us. After 40 ps from detonation, the distance traveled by 
the fragments is expected to be relatively linear with time. 
Considering the other trials, these parameters from the 


Table 4. Summary of fragment velocities as calculated from the Gurney Equation, and as determined by a least squares fit to a linear 


and exponential model (n=1), in ms'. 


Explosive Wall thickness Casing material Gurney (+12%) X-ray images 
linear fit exponential fit 

C4 6mm AISI 1026 1457 1275 1340 

TBX 1 3.8mm AISI 1026 1014 940 1025 

TBX 1 6mm AISI 1026 703 760 960 

TBX 1 6mm DCI 742 725 955 

TBX 1 6mm GCI 636 - (855)* 

TBX 3 6mm AISI 1026 1028 910 990 


*“ Estimated; several X-ray images of the GCI casings were not captured due to instrumentation problems. 
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Table 5. Parameters determined by least squares fit to exponen- 
tial model, with n= 1. For the model with n =2, the values for V, 
were nearly identical to the velocity determined by the linear 
model (after Table 4). There is no significant difference between 
the values of a and t between n= 1 and n= 2, except where noted. 


Explosive Wall Casing a T Vo 
thickness material Us Us mms” 
C4 6mm AISI 1026 9.21.2) :17.6 1.340 
TBX 1 3.8mm AISI 1026 19.5 26.8 1.025 
TBX 1 6mm AISI 1026 0° 60.2" 0.960 
TBX 1 6mm DCI 0° 49.5" 0.955 
TBX 1 6mm GCI - - (0.855)° 
TBX3 6mm AISI 1026 =.20.9 26.3 0.990 


* From the fit for n=2, reasonable values for a and t are approximately 
32 us. 

> From the fit for n=2, reasonable values for a and t are 26 us and 54 us 
respectively. 

° Estimated. 


exponential equations are consistent with the observation 
above that the cylinders are expanding at nearly constant 
velocity. 


4 Discussion 
4.1 Influence of Explosive Material 


The distance-time relationships for an AISI 1026 casing 
with a 6.0 mm wall casing as a function of explosive fill are 
shown in Figure 6. The C4 charge resulted in fragments 
having the highest velocity, due to the smaller size of the 
fragments, and the fact that C4 also generated the most 
product gas. Despite its run distance, the maximum frag- 
ment velocity for the TBX 3 series could be estimated from 
the X-ray images, as the casing expansion was delayed at the 
initiation end only. The other end received the full pressure 
loading of the detonation and appeared to expand and 
fracture in a manner consistent with the other series. The 
inclusion of RDX does improve the fragment velocity over 
those velocities seen in fragments generated by TBX 1. The 
larger time constants for the TBX fills, as presented in 
Table 5, could be attributed to the secondary combustion of 
the added metal. While some of the metal powder surely 
reacted within the casing, most of this reaction took place 
after the fragments had been propelled beyond this 
influence [2]. 

The experimentally determined velocities of the frag- 
ments for the C4 and TBX 3 fills are both about 10% lower 
than those calculated by the Gurney equation. The velocity 
of the fragments from C4 fills is key, as this provides a 
benchmark for these trials. The lower observed velocity 
from C4 is attributed to end effects, as the length to diameter 
ratio is about 2.5. The detonation and fragmentation 
behavior from C4, as seen in the X-ray images, is compa- 
rable to those of ideal (length to diameter ratio of greater 
than 6) casings. 

Although TBX 1 propelled fragments much slower than 
the other two explosives, the fragments were still propelled 
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Figure 6. Distance versus time for 6 mm thick AISI 1026 casings 
as a function of explosive fill. The lines are least squares fits. 


faster than predicted by the Gurney equation. This may be a 
result of metal igniting in the TBX 1 product gases prior to 
the rupture of the casing, adding enough heat to increase the 
gas energy to propel the fragments. As the Gurney equation 
assumes that the fraction of energy propelling the fragments 
of any charge is roughly the same, this result may indicate 
that a larger fraction of energy is available in liquid 
propellant/metal formulations. If this hypothesis were 
true, one would expect, then, to see an increase in the 
fragment velocities from TBX 3. There are three facts to 
consider. First, both C4 and TBX 3 contain RDX, where 
TBX 1 does not. Secondly, the amount of metal in TBX 3 is 
about half that of TBX 1. Finally, the casing ruptured sooner 
for TBX3 than for TBX 1. These factors indicate that 
TBX 3 behaves more like C4 than TBX 1, and that the 
casing ruptured before any energy could be added to the 
fragments from the ignition of the metal powder. 


4.2 Influence of Casing Thickness 


Figure 7 presents the velocity profiles for TBX 1 filled 
AISI 1026 steel casings of 3.8 mm and 6 mm wall thickness. 
The fragments from the thinner casing were observed to 
travel at roughly 940 m/s, which is 7% lower than the 
velocity calculated by the Gurney equation, and is a much 
higher velocity than for the 6 mm casing. These fragments 
had about the same velocity as those for the TBX 3 fill in a 
6mm thickwalled cylinder. It is noteworthy that TBX 1 
contains more metal (in this case magnesium) than TBX 3, 
which contains aluminum. The magnesium has a relatively 
higher affinity for oxygen than aluminum, forms a non- 
protective oxide, and has a higher vapor pressure. In 
comparison to the 6 mm thick casing with a TBX 1 fill, the 
thinner wall results in a faster expansion rate and lower 
confinement time. These points further support the hypoth- 
esis that a higher fraction of energy is available for 
propelling fragments. This fraction of energy is dependent 
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Figure 7. Distance versus time for AISI 1026 casings of two 
different wall thicknesses with TBX 1 fill. The lines are least 
squares fits. 


upon the confinement of the detonation products and the 
ignition delay of the metal powders. 


4.3 Influence of Casing Material 


The influence of the casing material can be seen in 
Figure 8 for TBX 1 fills. The velocity of the fragments from a 
DCI case is calculated to be higher than that for the AISI 
1026 steel. The experimental results show the opposite, but 
not to a great extent, as they are well within the error of 12%. 
Although in apparent agreement in Table 4, the results for 
GCT are less conclusive than those from DCI and AISI 1026 
steel, due to the shortage of points at the longer time delays. 
The effect of the casing material is not as great as the effect 
of the larger proportion of energy available to propel the 
fragments. 
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Figure 8. Distance versus time for various casing materials, all of 
6 mm thickness with TBX 1 fill. No lines are drawn to maintain 
clarity. 
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5 Conclusions 


The X-ray images captured the fracture behavior of the 
casings as a function of fill and material characteristics. The 
casings fragmented as expected. The X-ray images also 
provided information on the run-up distance of the explo- 
sive fills used here. The run distance for the TBX 3 
formulation, containing liquid monopropellant, aluminum 
powder, and RDX, is about 20 mm. The run distance for the 
TBX 1 formulation is less than 20mm, such that no 
indications of asymmetric expansion are observed. 

The Gurney equation assumes that the fraction of energy 
propelling the fragments of any charge is roughly the same. 
Results in this work are consistent with the conclusion that a 
larger fraction of energy is available in TBX (liquid propellant/ 
metal) formulations. Furthermore, this fraction of energy is 
dependent upon the confinement of the detonation products 
and the ignition delay of the metal powders. These two factors 
have a greater influence on the fragment velocities than do 
material characteristics of the casing. 
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Figure 9. A comparison of Gurney velocity to detonation 
velocity for various explosives, raw data taken from Cooper [5]. 


Appendix 
Calculating the Gumey Velocity 


Since the values of the Gurney constant for many 
explosives, including the TBXs in this research, are not 
available in the literature, it is important to establish a 
consistent method for their prediction. There is a roughly 
linear relationship between an explosive’s Gurney constant 
(V/2E) and VOD [5]. This is established as: 
VOD = 2.91V2E (Al) 

Figure 9 is based on the same data used by Cooper; 
however, the approximation is improved if the trend line is 
not forced through zero. This fit will not predict the Gurney 


constant as well as Cooper’s relationship for explosives with 
a very low VOD (<2000 ms‘). 
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Figure 10. A comparison of prompt to terminal values of the 
Gurney constant for various explosives, raw data taken from 
Dobratz [6]. 


Calculating the Prompt Gurney Velocity 


The Gurney constant is dependent on both the explosive 
and casing material. When brittle casings are used, the 
subscript p denotes the typically prompt casing fracture. If 
no subscript is shown, or the subscript is ¢ (for terminal), the 
constant denoted as V2E applies to casings that are ductile 
and fail primarily due to shear. Figure 10 compares the 
material dependency of the Gurney constant. This figure is 
based on work by Dobratz [6] that compares the Gurney 
constant for explosives in ductile casings and brittle casings. 
The linear relationship is used to approximate ,/2£E, for 
TBX 1 in a GCI casing. 
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